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1 Introduction 

In its most general sense, a (smooth) dynamical system, from the geometrical point 
of view, is simply a vector field on some manifold. But many dynamical systems of 
interest in physics and engineering apphcations are more speciahzed than that: they are 
of second-order type. By a dynamical system, or vector field, of second-order type, or a 
second-order differential equation field, we mean a vector field F on the tangent bundle 
r : T{Q) — )• Q of some configuration manifold Q with the property that r^yV = y for all 
y € T{Q), so that in terms of coordinates {x^,y^) where the are the canonical fibre 
coordinates corresponding to coordinates on Q, 

^ dx^ ~'~ dy^ 

It is of interest therefore, to find criteria for determining whether a given dynamical 
system, which may be represented in some arbitrary coordinates, is actually of second- 
order type, in that coordinates may be found with respect to which it takes the form 
above. This is a problem which has both a local aspect, just described, and a global one. 
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which includes the question of whether the manifold on which the dynamical system 
resides is in fact a tangent bundle. 

In the recent article [11], Ricardo and Respondek deal with a version of the problem in 
the context of control theory. Assume M to be an even dimensional manifold, not known 
to be the tangent manifold of some configuration manifold. Under which conditions does 
a coordinate change exist such that a given nonlinear control system 

= F"{z) +UrG^{z), a = l,...,2n, 

is transformed into the form of a so-called 'mechanical control system', meaning a dy- 
namical system of the form 

r = mQH'f + Pbi<l)Q' + Q^il) + urg^riz), a = 1, . . . , n? 

The solution of the problem in |11] is cast in terms of a certain vector space V, which is 
a subspace of the infinite dimensional vector space of vector fields on M, with dimension 
exactly the half of the dimension of the manifold, and which, among other properties, 
contains the control forces Gr and satisfies [V, V] = 0. 

In this paper, we will address a more general problem. First of all, we will not assume 
that the dimension of the manifold is even. This is mainly motivated by the observation 
that even for a system of time-dependent second-order differential equations the manifold 
on which the dynamics is described is odd-dimensional: it is the first jet manifold of a 
spacetime manifold, or event space, fibred over the real numbers (see e.g. [3l[7]). Further, 
we will allow some of the new coordinates to play simply the role of parameters. That 
is to say, we will not require the number n of second-order equations to be exactly the 
half of the dimension of the manifold M. 

Next to extending the results of [llj to a broader class of manifolds, we will also make 
some conceptual modifications. In a nutshell, the results of [TT] claim that if V (which is 
a vector space constructed from the given control forces Gr and from F) satisfies certain 
conditions the vector field F = F"d/dz°^ (the so-called drift vector field) transforms 
into an appropriate coordinate form, and, as a side-effect, so does also the controlled 
dynamical field F + UrGr- We will take the space V to be the primary given object of our 
study, and ignore that it was constructed from some given control forces. Consequently, 
we shift the attention to specific coordinate expressions for F, and leave the control 
system given by the vector field F + Uj-Gr out of the picture all together. A second 
deviation is that for us V will not be a vector space of vector fields, but rather the 
distribution it generates. 

In Section 2 we investigate under what conditions a given vector field F can be trans- 
formed into the coordinate expression of a second-order differential equation field with 
possible parameters, in the presence of an arbitrary involutive distribution V (of arbi- 
trary dimension). Our framework has the advantage that it leaves open the possibility 
that the transformed dynamics become either autonomous or time-dependent. We show 
in Section 3 how one can associate various connections to F, and we argue that these 
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connections provide a coordinate-independent method to express that the dynamics of 
F is of quadratic type (or of mechanical type, in the sense as above) in a yet unknown 
set of coordinates. 

Working with a distribution V has the further advantage that it brings an associated 
almost tangent structure (and almost jet structure) to the foreground. These geometric 
structures find their equivalence in standard tangent bundle and jet bundle geometry, 
but they went unnoticed in [TT]. Based on results in [H [SJ [TU] we further address in 
Section 4 the global issues that arise in this context, such as e.g. the affine fibre bundle 
structure of M and the relation of F to second-order differential equation fields on a 
certain tangent or jet manifold. 

In the last section we illustrate the theory in the context of a Lagrangian system with an 
Abelian symmetry group, where a second-order differential equation field with multiple 
parameters naturally shows up. 

2 Local coordinate transformations 

Let M be a manifold of dimension m and V an involutive distribution on M of dimension 
n such that 2n < m. 

We must be a little careful here about the meaning of the term distribution and related 
terms. A distribution on M is of course a choice of subspace of Tz{M) at each z G M, of 
constant dimension, depending smoothly on z in the sense that it admits local smooth 
bases. There is a related, but distinct, concept, which we may call a vector field system. 
A vector field system <S on M is a collection of (smooth) vector fields on M which is a 
C°°(M) submodule of X{M), the module of vector fields on M. For each z € M, we 
denote by dimz{S) the dimension of the subspace of Tz{M) spanned by the values at z 
of the vector fields in S. Now dimz{S) need not be constant. However, if vector fields 
Xa are linearly independent at z they are linearly independent in a neighbourhood of 
z, which means that dim^/(5) > dimz{S) for all z' in a neighbourhood of z. Moreover, 
dim^(5) has a maximal value on M, which we call the maximal dimension of S, and the 
set of points at which the maximal dimension is attained is an open subset of M. 

A vector field system 5 on M restricts to a vector field system S\u on any open subset 
U of M, considered as a submanifold. On the other hand, any vector field in Sjj may be 
extended to a vector field in S by multiplying it by a bump function whose support is 
contained in U (taking advantage of the fact that we are working in the C°° category). 
So it is permissible to discuss local aspects of vector field systems in coordinates. 

An alternative definition of the term distribution is that a distribution is a vector field 
system S for which dim2(5) is constant, or for which the maximal dimension is attained 
everywhere. 

We assume that V, mentioned in the opening sentence of the section, is a distribution in 
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the strict sense. Now suppose that we have a vector field F on M not belonging to V. 
We denote by V + [F, V] the collection of vector fields on M which may be written in the 
form Vi + [F, V2] with Vi, V2 G V. This is a vector field system, essentially because for 
any / G C~(M) and V eV, f[F, V] = [F, fV] (mod V). We will be concerned with this 
vector field system, for a given involutive distribution V and vector field F, throughout 
this paper. 

Proposition 1. Suppose that the vector field F is such that [F,V] n V = {0}, that is, if 
V (zV and [F, V] then V = 0. Then the maximal dim^ension 0/ V + [F, V] is 2n, and 
the open subset of M on which it is attained is dense, that is, its closure is M. 



Proof. We denote the vector field system V + [F, V] by W for convenience. Clearly the 
maximal dimension of W is at most 2n, and the set of points z where dim^Cl/V) = 2n is 
open, though it may be empty. Suppose that z is a point of M with dim2(W) < 2n. We 
show that there can be no open neighbourhood of z such that dim2;/(W) < 2n for all z' 
in the neighbourhood. 

Since V is involutive there is a coordinate neighbourhood U of z and coordinates (g", 
with a = 1,2, . . . ,m — n and i = 1, 2, . . . , n such that the coordinate fields d/dy'^ span 
V\u- We may set 

d d 



for some smooth functions f^, on U, so that 



F\u 



d_ 

' dy^ 



Qyi Qqa 

and for any V G Vu, with V = V^d/dy^, 



(mod V), 



Now since m — n> n the rank of the matrix {df^ jdy^) at any point is at most n. If it 
is n at z then 

V\z)^{z) = Q =^ V\z) = Q, 

the vector fields [F\u,d/dy'^] are linearly independent at z, and dim^CVV) = 2n. So if 
dim2;/(>V) < 2n for all z' in a neighbourhood of z, which we can take to be a coordinate 

neighbourhood as above, then the rank of the matrix {df'^/dy^){z') is less than n, and 
we can find functions on a neighbourhood U of z, not all vanishing, such that 

Then the vector field V = V^d/dy^ on U satisfies [-F|(7,y] G V. So (by multiplying by 
a suitable bump function) we can find a vector field V on M, not identically zero, with 
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V' G V, such that [F, V'] € V, which is a contradiction. So every neighbourhood of a 
point z where dim^CW) < 2n must contain a point z' where dim^/(W) = 2n. Thus the 
set of points z where dim^CW) = 2n is nonempty and open, and its closure is M; that 
is, the set of points where dim^Cl/V) = 2n is an open dense subset of M. □ 



From now on we will assume that the vector field -F does indeed satisfy the condition of 
Proposition [H and we will restrict our attention to the open subset where dim^CW) = 2n, 
that is to say, we will effectively assume that W is a distribution. We will make the further 
assumption that the distribution W is also involutive. 

We will work locally for the rest of this section, and drop explicit notational reference to 
the neighbourhood on which we are working. 

In the proof of Proposition [T] we showed that if {Vi} is a local basis of V consisting of 
coordinate fields d/dy^ of a local coordinate system and we set Wi = [F,Vj\ then 

{Vi, Wj} is a local basis for W. Indeed, this will be true for any local basis {V^} of V. 
If we change basis to = A-'^Vj (where the A-'^ are locally defined smooth functions and 
{AD is nonsingular) then Wj changes to Wi = AjWj + F{Al)Vj. If the basis {Vi} is 
such that [Vi,V^] = 0, so that the Vi are coordinate fields, the necessary and sufficient 
condition for {V} also to satisfy [Vi,Vj] = is that A{Vi{Af) = A^jVi{Af). 

Our next aim is to show that, under the assumptions stated earlier, one can choose a 
commuting basis {Vi} for V such that it and the corresponding Wi satisfy [y^jH^j] € V. 
For any basis {Vi} we can write 

[V,Wj] = a'ljVk + P^jWk. 
Notice that if [V^, V^] = then both coefficients are symmetric in their lower indices: 

= [F, [V, Vj]] = m, Vj] + [Vi, W,] = [V, Wj] - [Vj,Wi]. 
If we change the basis of V to = A-^- Vj we have 

[Vi, Wj] = [A\Vk,A]Wi + F{A])Vi] 

= [A'iViiA]) + A\A^pfjWk (mod V). 

So to make [Vi, Wj] G V we want to choose such that Vi{A'^) + AfPj'^ = 0. Note that 
since is symmetric, we will then have A[Vi{Aj) = AjVi{A^), and if the Vi pairwise 
commute then the Vi will also pairwise commute. 



The equations 



dy 



^ Plm 



are linear first-order partial differential equations for the unknowns A^-, and admit solu- 
tions if and only if their integrability conditions, which are 

Qyi Qyj ~^ PimPjk PjmPik ^) 
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are satisfied. Now 
and 

[Vi, [Vj,Wk]] = mp'j,) + PLPTkWl (mod V). 

So it follows from the identity — [V}, [V^, PF^]] = that the integrability 

conditions are indeed satisfied. If we take a solution Aj for which the matrix (Aj) is 
nonsingular on a local cross-section of the V foliation, for example by taking Aj = 5j 
there, then (A^j) will be nonsingular on an open subset containing the cross-section. We 
have shown the following. 

Proposition 2. // both V and W are involutive, there is a commuting basis {Vi} of V 
such that for all i,j, [Vi,Wj] G V (where Wi = [F,Vi]). 

The remaining freedom in the choice of commuting basis (such that [Vi,VFj] G V still 
holds) is to take A'- to satisfy T4(^|) = 0. 

The condition [Vi, Wj] G V says that Wj is invariant under the action of V, modulo V. 

Let us take a coordinate neighbourhood U in M, with coordinates (g",y*) such that 
Vi = d/dy''; we may suppose without essential loss of generality that U is the image 
of a product of open subsets O C M'""" and P C M", where £ P. Then = is 
a submanifold of U of codimension n, say N, and U is fibered over N with fibres the 
integral submanifolds of V. Denote hy n : U —?■ N the corresponding projection. Then 
the restriction of Wj to U is projcctablc to N: that is to say, there is a well-defined 
vector field Wj on N which is 7r-related to Wj. More generally, a vector field X gW, say 
X = X'^Wi (mod V), is projectable if, and only if, the coefficients X'^ satisfy V{X^) = 
for all y G V, or indeed if Vj{X^) = 0. Let us denote by VV the distribution on spanned 
by the Wi, in other words, the distribution consisting of the projections of projectable 
vector fields in W. Then VV is involutive, since it 7r-related to the involutive distribution 
W. We may therefore choose coordinates (t^, x*) on TV, where p = 1, 2, . . . , m — 2n, such 
that the integral submanifolds of VV are given by t^ = constant. Then with respect to 
the coordinates {tP,x^,y^) on U we have 

^^ = ^' W^. = Wf(x)A („,odV), 

where the coefficients Wf are everywhere the components of a nonsingular matrix. We 
still have at our disposal the freedom to change the original basis to = A'^{f,x)Vj. If 
we do so with A^Wj, = then 

If we make a further change of coordinates to 

iP = tP, x' = -x\ f = Wi{x)y^, 
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then 

We have proved the following result. 

Proposition 3. With V, W, Vi, Wi as above, we can find local coordinates {tP,x^,y^] 
on M, p = 1,2, . . . ,m — 2n, i = 1,2, . . . ,n such that 

We move on now to investigate the form of F. There are two cases to consider, depending 
on whether F does or does not belong to W. 

Theorem 1. Assume both V and W are involutive. 



1. Suppose that F G W, and assume that the set N C M = {z & M : F{z) € V} is 
nonempty. Then we may choose coordinates with respect to which 

d d 

F = y'^+r{t,x,y)^; 

that is to say, F takes the form of a second- order differential equation field in terms 
of the coordinates with the merely behaving as parameters. 

2. Suppose that F is everywhere independent of W (so that in particular M > 2n) 
and that [F, W] C W. Then we may choose coordinates with respect to which 

that is to say, F takes the form of a time- dependent second- order differential equa- 
tion field in terms of the coordinates {t^,x^,y'^), with the t^ with p > 1 merely 
behaving as parameters. 



Proof. 1. For the first case, set F = a^Vi + UWi with respect to a frame with [V^,, Vj] = 
and [Vi,Wj] e V. Then 

Wi = [F, Vi] = -V{a^)Vj + V[Wj, - V,{W)Wj, 

so we must have Vi{bP) = Now is the zero level set of (6*), and the rank of the 
Jacobian of the map M R" : z i— )• (^^(-z)) is n, or in other words the 1-forms db^ are 
independent, since {Vi{V)) is nonsingular. So is an immersed submanifold of M of 
codimension n, and V is transverse to it. We may choose coordinates as in 

Proposition [31 such that Vi = d/dy^, N is given by = 0, are coordinates on A^, 

and 

Wi = -^ (modV). 
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With respect to such coordinates set 

F = f\t,x,y)— + F'{t,x,y) — . 

Then from its definition 

d 

and therefore 

Taking into account the fact that P{t,x,0) = we have f^{t,x,y) = y*. 

2. For the second case, take coordinates as in Proposition O and suppose that 

F = ^P{t,x,y)^ (modW). 

From the assumption that [F, W] C W it follows that in fact (fP depends only on the 
By assumption the ip^ cannot vanish simultaneously, and so by a transformation of the 
coordinates we may take 

Arguing as above we see that 

dy' 

We may only conclude now that x, y) = + x). However, a further coordinate 
transformation y^ y^ + k^{t,x), with and unchanged, leads to 

d d d d 

—— 1-^ —— (mod V), 7-^ i-> 7-^ (mod V), 

dtp dtp ^ ' dx' dx^ ^ ' 

with d/dy^ unchanged, and so leads to the required form for F. □ 

The remaining freedom in transforming the coordinates (x*,y') in the first case, so as to 
preserve the form of -F, is 

3x^ 

x' = x\x), ff = —yK 

That is to say, the transform like canonical fibre coordinates on a tangent bundle. 
From this point of view it is natural to think of the coordinates in use before the final 
transformation leading to Proposition [3] as quasi- velocities. 

The remaining freedom in the second case is 

Here the transform like the jet coordinates of the 1-jet bundle of a manifold fibred 
over R. 
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3 Induced connections 



The coefficients (3^^ used in tlie proof of Proposition [2] liave the appearance of the compo- 
nents of a connection, and the integrability conditions quoted in the proof have the form 
of the vanishing of the curvature of this connection. We begin this section by explaining 
in what sense the f3^j are indeed the components of a flat symmetric connection. 

Let V be an involutive distribution on any manifold M. For any vector field X on M 
denote hy X + V the equivalence class of X modulo V, that is, the collection of vector 
fields differing from X by an element of V. For any 1/ € V, set 

DviX + V) = [V,X] + V. 

This is a well-defined operator on equivalence classes, which is M-linear in both argu- 
ments, and for / G C°°{M) satisfies 

Dfv{X + V) = fDviX + V), DvifX + V) = fDv{X + V) + V{f){X + V). 

That is to say, D has connection-like properties. By the Jacobi identity, for any Vi,V2 G V 

Dv.Dv, (X + V) - Dv.Dv, [X + V]- D^y^y^^ {X + V) = 0; 

that is to say, if D were a connection it would have zero curvature. 

More particularly, let W be another involutive distribution on M, with V C W; then 
we may restrict X in the construction above to lie in W. The same conclusions hold, 
mutatis mutandis. We may think of V and W as vector sub-bundles of T(M), and vector 
fields in the distributions as sections of the corresponding bundles V M, W M. If 
Vl^ is a section of W ^ M then we may think of TV -|- V as a section of the vector bundle 
W/V — 7- M. Then (using the terminology of Lie algebroid theory) Z) is a V-connection 
on W/V. 

Now take dim W = 2n, dim V = n, and suppose there is a type (1, 1) tensor field S" on W 
(that is, a section of the bundle W^W* — > M) with the algebraic properties of an almost 
tangent structure (so that imS* = kerS), with kernel V. Then S defines an isomorphism 
between sections of W/V and sections of V. So we may define a V-connection on V, say 
V,by 

Va = S{Dv, {W + V)) for anyW eW such that S{W) = V2. 

That is, VvS{W) = S{Dv{W + V)) = S{[V, W]). This is well-defined as a V-connection, 
and has vanishing curvature. Since V is a V-connection on V, it makes sense to talk 
about its torsion. But for any Wi,W2 S W, 

= S[SiWi), W2] - S[SiW2),Wi] - [SiWi),SiW2)] 
= -{[S{Wi), S{W2)] - S[S{Wi),W2] - S[Wi,S{W2)]) 
= -Ns{Wi,W2). 
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That is to say, the torsion vanishes if and only if the formal Nijenhuis torsion Ns of S 
(a type (2, 1) W-tensor) vanishes. 

If S has vanishing Nijenhuis torsion in this sense, and we restrict attention to any leaf 
of the involutive distribution V, we obtain a flat symmetric connection there. 

We now show how to construct such a W-tensor S in the case of interest. 

Proposition 4. Assume both V and W are involutive. There is a unique type (1, 1) 
W-tensor field S for which 

S{V) = and S{[F, V]) = -V, V e V; 

it satisfies ker S = imS = V and Ns = 0. 

Proof. Let {Vi} be a basis of V, and set Wi = [F, Vj]: then {Vi, Wi} is a basis for W. So it 
is enough to know how S acts on elements of the form V and [F, V]. The definition above 
is consistent: if V E V, then also fV G V with / a function on M, and S{[F,fV]) = 
S{F{f)V) + S{f[F,V]) = fS{[F,V]) = -fV. We have = Q, kev S = imS = V. The 
formal Nijenhuis torsion Ns obviously vanishes for two elements in V. Moreover, for any 
Vi,V2eV, 

Ns{Vi,[F,V2]) = -S[Vi,V2] = 

because of the assumed integrability of V. Likewise, by making use of the Jacobi identity 
(and because [Vi, V2] G V), 

Ns{[F, Vi], {F, V2]) = [Vi, 1/2] + S[Vu [F, V2]] + S[{F, Vi], V2] 

= [VuV2] + S[F,[VuV2]]=0. □ 

In the case of interest, where S is as defined above, the P^j are the connection coefficients 
of this connection with respect to the basis {Vi}. 

Suppose we have a further distribution T-L on M, of dimension n, contained in W, and 
everywhere transverse to V; in other words a complement to V in W. We call such a 
distribution horizontal. Then the restriction of 5" to 7^ is a C°°(-M)-isomorphism "H ^ V. 
For any V eV, denote by the unique element of Ti such that S{V^) = V. 

We can extend the V-connection V on V to a W-connection on V as follows: for any 
W and any F G V set 

VwV = Pv{[Pn{W), V]) + S{[Pv{W), F"]), 

where and Py ai'e the projectors on % and V, respectively. The right-hand side 
belongs to V and depends R-linearly on the arguments. For / G C°°(M) 

VfwV = PviifPniW), V]) + S{[fPv{W), n) 

= f^wV - V{f)Pv{Pn{W)) - F"(/)5(Pv W) 



10 



while 



VwifV) = Pv{[Pn{W), fV]) + S{[Pv{W), /n) 
= fVwV + PH{W){f)Pv{V) + Pv 
= f^wV + {Pn{W) + PviW))if)V 
= fVwV + W{f)V. 

So V is a covariant derivative. 

If G V, say W = Vi, then the new definition gives VV1V2 = 5([Vi,y2"]). According 
to the old definition, Vv^V2 = S{[Vi,W]) for any W such that 5(1^) = V2. But 1^ = 
is such that S{W) = V2; so the two definitions agree in this case. On the other 
hand, suppose that W £ H and that W is projectable in the sense that [W, V] C V 
(the horizontal projection of any projectable vector field is projectable, and the Wi are 
projectable as we pointed out before). Then = [W,V]. 

Assuming as before that [F, W] C W (which is automatically the case F £ W, and is 
an assumption in Part 2 of Theorem [1] if not), it is possible to define a Lie derivative by 
F of W-tensors: for example, in the case of a type (1, 1) W-tensor T as the commutator 
of operators ad F and T: 

{CfT){W) = [F,T{W)]-T[F,W]. 
Proposition 5. // [i^, W] C W, the vector field F defines a complement TiofVin W. 

Proof. We show now that, with the above definition, CpS defines two projection opera- 
tors on W. We first show that {CpS)^ = id. We have, for V £V, 

{CfS){V) = [F,S{V)]-S[F,V] = V 

and thus {CpSfiV) = V. Also, 

{CfS){[F, V]) = [F, S[F, V]] - S[F, [F, V]] = -[F, V] - S[F, [F, V]]. 

Since S[F, [F, V]] € V, we have {{CfS){S[F, [F, V]]) = S[F, [F, V]], and therefore 

iCFSfi[F,V]) = -iCpSmV]) - S[F, [F,V]] = [F,V]. 

The conclusion is that P-^ = ^{id —CpS) and Py = ^{id+CpS) are complementary 
projection operators, with e.g. Pv{V) = V and P^(y) = 0; Ti = imP^ is therefore a 
complement to V in W. □ 

We will use this complement from now on. 

In the case where M is a tangent manifold T(Q), we can take V to be the canonical vertical 
distribution, and in particular W = X{T{Q)). The connection with covariant derivative 
V is then the Berwald connection associated to a system of autonomous second-order 
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differential equations, see e.g. [2] (taking into account the fact that the current connection 
is expressed in terms of vertical vector fields rather than vector fields along the tangent 
bundle projection). A similar construction exists for the case where M is the first jet 
manifold of a bundle E' ^ R, and where the second-order dynamics are time-dependent, 
see e.g. [3 E]- 

The Berwald connection can be used to describe special classes of second-order differential 
equation fields, such as the ones of quadratic type we had encountered in the introduction. 

The 'mixed curvature' of the V-connection V is the (1,2) V-tensor field 6 given by 

Proposition 6. Let V and W both he involutive and [F, W] C W. The necessary and suf- 
ficient condition for the existence of coordinates in which F takes the form of a quadratic 
second- order differential equation field is that 6 = 0. 

Proof. Let V denote the set of y G V for which the corresponding is projectable, i.e. 
V satisfies [1^",V2] G V for all V2 G V. This set defines a module over the projectable 
functions on M (those functions / for which Vi{f) = 0, for all Vi G V). Alternatively, 
^2 G V if and only if Vyi^a = 0, for ah Vi G V. 

Let Vi,V2£V. Then 

Vv,VyuV2 = -9{Vi,V3)V2 - V[v'Hy3]y2 + VyuVv,V2 = -9{Vi,V3)V2, 

meaning that VyiiV2 is again projectable if and only if 9 = 0. 
Let 9 = 0. If we set, for Fi, V2 G V, 

Dv,V2 = VyH^a, 

one easily verifies that the operator D satisfies connection-like properties with respect to 
the multiplication of elements of the module V with projectable functions /: 

Dfv,V2 = fDv,V2 and DvJV2 = fDv,V2 + Vf{f)V2. 
Remark that V^{f) is again projectable, since for any V2 G V, V2V^{f) = V^V2{f) + 

[V2,vm) = o- 

In the coordinates as defined in Theorem [1] (regardless of whether F lies in W or not) 
the connection coefficients are given by 

d d BP dT^- 

It is clear from this expression that T^j = Tj-, or, equivalently, that the connection D is 
symmetric, in the sense that the torsion 

Dv,V2-Dv,Vi-S[V^'',V2^] 
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vanishes. 

We can therefore conclude that the functions T^j are projectable if and only if = 0. 
For that to be the case, f'' must be of the form = Vf^y^y^ + P^y^ + , for some 
projectable functions Pj'{t,x) and Q^{t,x). □ 

The advantage of the current description is that the criterion = can be verified in 
any given set of coordinates on M. 



4 Global properties 

In this section we will address the global bundle structure of a manifold M in the context 
of a given involutive distribution and a vector field F, assuming from the start that the 
set where the dimension of V + [-F, V] is maximal is the whole of M, or in other words 
that V + [F, V] is actually a distribution. 

Let M be an m-dimensional manifold with an involutive (and thus integrable) n-dimensional 
distribution V. The foliation of the distribution defines an equivalence relation on M 
by declaring two elements of M to be equivalent if they lie in the same leaf of V. The 
quotient of M by means of this equivalence relation, say Q, will have the structure of a 
differentiable manifold if for every leaf one can find a smooth embedded local submani- 
fold N through a point of the leaf, of dimension m — n, which has the property that each 
other leaf it meets is intersected in only one point. Then vri : M — >■ Q defines a fibration, 
for which the fibres have dimension n, and for which the distribution V coincides with 
the tangents to the fibres. 

We will assume that this condition is satisfied. 

In the case of interest, Q comes also equipped with an integrable distribution. Indeed, 
a projectable vector field on M is vri-related to a vector field on Q. Those projectable 
vector fields that happen to lie in W define therefore a distribution on Q, say W, which 
is involutive by construction. As above, the corresponding equivalence relation therefore 
defines a new quotient, T, again assumed to be a manifold. We will denote the corre- 
sponding fibration by 7r2 : Q ^ T. Alternatively, we could have defined a fibration by 
quotienting out the distribution W from the beginning. This structure will coincide with 
the composition projection 7r2 o vri : M ^ T. 

An almost tangent structure on an even dimensional manifold is a (l,l)-tensor field S 
on that manifold, for which the kernel of S at each point coincides with its image. The 
almost tangent structure is said to be integrable if its Nijenhuis torsion vanishes. If that 
is the case then the kernel of S is an involutive distribution. We now recall a result from 
[4]. Suppose that the kernel of an integrable almost tangent structure on a manifold 
defines a fibration over some base manifold (as above, by taking the quotient of that 
distribution). Suppose that each fibre is connected and simply connected, and that there 
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exists a flat connection on each fibre, for which the fibre is geodesically complete. Then 
the manifold is the total manifold of an affine bundle modelled over the tangent bundle 
of the base manifold. 

This theorem can be applied to the current setting, if we take a particular leaf of W 
to be the even dimensional manifold of interest. The projection tti will project this leaf 
Lyy of W onto a corresponding leaf Lyy of VV. Therefore we may consider the fibration 
given by the restriction vri|i^ : Lw — > ^Vv- Vectors that are tangent to its fibres can 
be identified with vectors in V, and the fibres themselves can be identified with leaves 
of V. We have defined a W-tensor field S" on W in Proposition HI It restricts naturally 
to an almost tangent structure on Lyy (i.e. S'^ = 0, Ns = 0, and kerS" = im5 = Vl^). 
Moreover, we have seen that the restriction of the V-connection Vv {V G V) to a leaf of 
V gives a flat connection on that leaf. We can conclude therefore: 

Theorem 2. Suppose that each leaf of V is connected and simply connected, and assume 
that each leaf of V is geodesically complete with respect to the restriction of Wy to that 
leaf. Then, for each Ly^ ofW, ttiIl^ : Lw Ly^ is an affine bundle, modelled over the 
tangent bundle T{Lyy) — )• Lyy. Suppose further that the set N = {z G Lyy : F{z) € V} is 
a global cross-section of ttiIl^: then Ly\; may be identified with T{Ly^) and N with the 
zero section. 

Corollary 1. In case F S W, and under the assumptions of the previous theorem, the 
restriction of F to a certain leaf Ly^; = T{Ly;) will be a second-order differential equation 
field on T{Ly^). 

Proof. This follows easily from the coordinate expression of F. Restricting F to a leaf 
is the same as fixing the parameters t^ to some constant values. □ 

For completeness, we mention that one may find an alternative formulation of the theo- 
rem of [1] in [S], where the global conditions on an (assumed given) symmetric connec- 
tion are replaced by global conditions on an (assumed given) vector field. In our current 
framework, the restriction of the vector field S{F) to L>v plays the role of that vector 
field. 

In case F does not belong to W, but leaves it invariant, F defines a vector field F on T. 
This vector field defines a 1-dimensional involutive distribution on T, leading as before 
to a fibration T — >■ Tq. If we assume that the vector field F is complete, an integral curve 
of F will define a 1-dimensional submanifold Ti of T. In turn, the preimage of Ti under 
is a collection E of leaves of VV lying over that integral curve. We can think of the 
restriction of tt2 to E as defining a fibration 7r2|_B : E ^ Ti. Let's denote its 1-jet bundle 
by J^{E) E. 

The distribution yVp = (-F)©W is also involutive. Its leaves are (2n-|-l)-dimensional 
manifolds that are projected by means of vri onto one of the above described manifolds 
corresponding to a certain integral curve (with image Ti) of F. The fibres of vril^^ : 
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~^ ^ again be identified with V. Recall that we had defined a symmetric flat 
connection on each leaf of V. 

By setting S{F) = we can extend 5 to a Wir-tensor field, which has the property that 
its restriction to a leaf of Wf satisfies S*^ = 0, Ns = and rank 5 = n. These properties 
are exactly those that define, in the terminology of |10) . an 'almost jet structure'. The 
above mentioned theorem in [3] has been generalized to (2n + l)-dimensional manifolds 
with almost jet structures in pXi], where a 1-jet bundle replaces the role played by a 
tangent bundle (see the Theorem on page 90 of [10]). We are in the situation that we 
can apply this theorem, since each Lyy^ is a (2n + l)-dimensional manifold with a 2n- 
dimensional distribution W|_Lyy^ which is completely integrable and which is such that 
imS" = Vjiyy^ C W[Lyy^. Wb 1018^ thcrcforc conclude: 

Theorem 3. Suppose that each leaf ofV is connected and simply connected, and assume 
that each leaf of V is geodesically complete with respect to the restriction of Vy to that 
leaf. In case F and [F,W] C W, each leaf L-^p ofWp is diffeomorphic to the 1-jet 
bundle J^{E) ofE^Ti. 

Corollary 2. Under the assumptions of the previous theorem, the restriction of F to 
a certain leaf L\^p = J^{E) will be a time- dependent second-order differential equation 
field on J^{E). 

Proof. This follows again from the coordinate expression of F. Restricting F to a leaf is 
the same as fixing all parameters t^ to some constant values, except for t^. □ 

The cases of most obvious interest are those in which the dimension of M is either 2n or 
2n + 1 (n being the dimension of V) . We end this section with a statement of our global 
results in these cases, in a form which collects together the assumptions we have made. 

Theorem 4. Let V be an involutive distribution of dimension n on a manifold M of 
dimension m, m = 2n or 2n + 1; and F a vector field such that VH [-F, V] = {0}. Assume 
that 

• M is fibred over a manifold Q and the leaves of V are the fibres of this fibration; 

• each leaf of V is connected and simply connected; 

• each leaf of V is geodesically complete with respect to the flat symmetric connection 
induced on it (as described in Section 3); 

• V + [-F, V] is a distribution (necessarily of dimension 2n). 
In the case m = 2n, assume further that 

• the set {z G M : F[z) G V} is a global cross-section of M ^ Q. 
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Then M may be identified with T{Q) and F with a second-order differential equation 
field on T{Q). 

In the case m = 2n + 1 assume further that 
• W = V + [-F, V] is involutive; 



• F is complete; 

• Q is fibred overW, and the leaves of W (the projection of W to Q) are the fibres 
of this fibration. 

Then M may be identified with J^{Q) and F with a time- dependent second-order differ- 
ential equation field on J^{Q). 

5 An illustrative example 

The cases in which dimM is 2n or 2n + 1 may be of most obvious interest, but they are 
by no means the only cases of interest, as we now show by an example. 

Let Q be the configuration space of a Lagrangian system with regular Lagrangian L 
and assume that L is invariant under the (free and proper) action of a symmetry Lie 
group G. In that case, the Euler-Lagrange field T G X(T(Q)) is G-invariant and it can 
be reduced to a vector field T on T(Q)/G. The corresponding equations for finding the 
integral curves of T are known in the literature as the 'Lagrange-Poincare equations', see 
e.g. [1]. The equations determining the reduced vector field F € X(T(Q)/G) can be cast 
in terms of the reduced Lagrangian / on T(Q)/G. 

We will follow here closely the description we have given in [8] . With the aid of a principal 
connection on Q Q/G one may decompose T{Q)/G into T{Q/G) ® {Q x g)/G, where 
the action of G on g is the adjoint action. In what follows {x^,v^,wP) are local coordinates 
on T{Q)/G, where the (x*) are coordinates on Q/G, and {v^,w^) are fibre coordinates 
corresponding to the decomposition. 

We will assume that the symmetry group is Abelian. This has the advantage that the 
adjoint action is trivial. The vector field T can then be determined from: 



F ^W, [F, W] C W; 
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Here Kfj^ are the components of the curvature of the principal connection with respect 
to an invariant basis. The coordinate expression of the reduced field is therefore of the 
form r = v^d/dx^ + T^d/dv^ + T^d/dw^ , where and are functions on T{Q)/G. 

In the assumption that the matrix {d'^l/dw^dw'^) is everywhere non-singular, and that 
the relation dl/dw^ = jip can therefore be rewritten in the form = pP{x,v,fj,), we 
can perform a coordinate transformation — )• (x* = x^,v^ = v^,^p = dl/dw^). 
The last equation is then simply T {fj,p) = 0, that is to say: the coordinates /ip can 
be regarded as parameters. In the new coordinates the reduced vector field becomes 



The first two equations determine a system of second-order ordinary differential equations 



7^^(x, v) = l{x, V, p{x, V, fi)) - HpfP{x, V, /i) 
these equations, restricted to fixed values for the fis, can equivalently be rewritten as 



This equation is known as Routh's (reduced) equation for an Abelian symmetry group, 
see e.g. [6]. 

We show that the situation described above is in agreement with the statements of 
Theorem[T]and Corollary[TJ Recall first the definition of the momentum map J : T{Q) 
g*, where (J(f),^) = {S,^L){v) (for each £ Q, S, £ X(Q) is the corresponding fundamental 
vector field). It is well-known that the map T(Q) — )• Q x Q*,{q,v) i-^ {q,J{v)) is G- 
equivariant, where the action of G on g* is the coadjoint action. Therefore, it reduces 
to a map J : T{Q)/G {Q x g*)/G. But for an Abelian group the adjoint action is 
trivial, so the coadjoint action is also trivial. It follows that J is invariant, and that the 
image of the reduced momentum J is Q/G x g*. In the current coordinates J is simply 
{x,v,n) ^ {x,fj.). 

Let M = T{Q)/G and F = T. The distribution V = kerTJ is clearly involutive. 
It has the commuting basis given by the vector fields d/dv^. It is easy to see that 
[F, V] n V = {0}. The distribution W = V + [F,V] is spanned by {d/dx\ d/dv'} and is 
involutive as well. It is the distribution formed by the vector fields on T(Q)/G which 
are tangent to the level sets of momentum. A leaf of W is thus a particular level set, 
fip = fip. The corresponding N in the statement of Theorem [1] can be identified with 
im J = Q/G X g*, and is non-empty. So Theorem [1] applies. 

The quotient space Q of Section H] can be identified with Q/G x g*. It is trivially fibred 
over T = g*; VV is the distribution formed by the (projections of the) d/dx^, and any 
leaf of VV can therefore be identified with Q/G. According to Theorem [2l therefore, the 
restriction of to a level set of momentum ji = fi^ is a second-order differential equation 
field on Q/G (a leaf of VV); it is of course the one which satisfies Routh's reduced equation 
with fip = fip. 



r = v'd/dx' + rd/dv' + Od/dnp. 



in the variables with parameters fip. By introducing Routh's (reduced) function 
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